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Abstract 

This  report  documents  an  analytical  model  to  assist  in  the  planning  of  the  Operational 
Test  and  Evaluation  (OT&E)  of  the  Marine  Corps'  prospective  Advanced  Amphibious 
Assault  Vehicle  (AAAV).  The  model's  emphasis  is  on  suitability  issues  such  as 
Operational  Availability  in  an  on-land  (after  ocean  transit)  mission  region.  The  model 
predicts  strong  sensitivity  to  the  form  of  an  assumed  distribution  of  times  to  vehicle 
breakdown,  and  hence  recommends  that  appropriate  test  data  be  obtained  to  reveal  that 
form  (the  mean  alone  is  inadequate).  Removal  of  design  faults  likely  to  cause  early 
failures  is  encouraged. 

The  AAAV  design  is  for  a  relatively  lightweight  but  technologically  advanced 
(mobile  and  lethal)  amphibious  vehicle  that  operates  in  platoons.  If  individual  vehicles 
break  down  they  must  often  be  transported  for  repair,  for  example,  to  their  point  of 
origin,  a  ship.  One  option  is  to  allow  another  platoon  member  to  tow;  another  is  to  assign 
an  auxiliary  vehicle  to  transport.  Such  duties  may  seriously  diminish  the  platoon 
productive  mission  availability  (the  towing/transport  agents  may  themselves  fail).  OT&E 
should  be  designed  to  fully  test  for  the  system-wide  effect  of  force  (platoon,  and  beyond) 
sustainability  requirements. 
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Executive  Summary 


This  report  documents  dynamical  and  probabilistic  models  to  assist  in  planning  the 
Operational  Testing  of  the  Marine  Corps’  Advanced  Amphibious  Assault  Vehicle  (AAAV).  The 
model  is  programmed  in  rapidly  executable  software  available  from  the  authors;  it  runs  on  a 
personal  computer  (PC  to  readers  of  this  report).  The  software  allows  an  analyst  to  quantitatively 
study  the  sensitivity  of  the  operational  availability  of  one  (or  more)  platoons,  each  of  12  vehicles. 
Among  sensitive  features  or  properties  are: 

•  The  statistical  nature  (inferred  distribution,  not  just  the  “mean”)  of  the  times  from  a  repair 
completion  to  subsequent  failure.  For  sensitivity  to  distribution  form  see  Sections  4  and  6. 

•  The  times  to  repair  failed  vehicles  (including  the  often  significant  time  and  resources 
required  to  transport  those  vehicles  to  the  repair  facilities). 

•  The  numbers,  types,  and  concepts  of  employment  (COE)  of  the  necessary  support 
vehicles  (e.g.,  landing  craft  and  helicopters). 

The  motivation  is  this:  modem  military  equipments  tend  to  be  designed  to  be  lighter,  more 
transportable,  and  quickly  mobile,  but  also  more  combat  effective  than  their  predecessors.  An 
example  is  the  Marine  Corps’  AAAV,  which  is  intended  to  replace  the  present  Amphibious 
Assault  Vehicle  (AAV).  Although  it  is  smaller  than  the  AAV,  the  new  AAAV  is  faster  and  more 
technologically  advanced.  It  possesses  more  sophisticated  sensors  and  more  lethal  weaponry  and 
is  expected  to  be  more  combat  effective  than  the  present  AAV. 

Such  technological  improvements  are,  however,  likely  to  be  characterized,  at  least  initially, 
by  increased  fragility,  and  certainly  less  self-contained  on-board  sustainment  (repair)  capability. 
For  the  AAAV,  the  advanced  technology  that  results  in  greater  speed  and  lethality  comes  at  the 
expense  of  external  support  requirements  in  case  of  onboard  equipment  failures  (e.g.,  of  sensors 
and  communications  gear),  and  particularly  in  case  of  mechanical  breakdowns.  In  general. 
Operational  Testing  of  new,  sophisticated,  potentially  highly  capable  and  effective,  but 
potentially  failure-prone  assets  should  be  done  in  the  context  of  their  entire  essential  support 
system.  Such  tests,  if  done  entirely  in  the  field,  are  costly  and  time-consuming  at  best. 
Consequently,  a  trustworthy,  model-based  preview  of  the  entire  system  in  operational  action 
should  be  of  great  suggestive  value.  Focused,  smaller-scale  tests  along  with  appropriate  data 
acquisition  and  analysis,  can  be  used  iteratively  to  parameterize  overall  models;  such  Model-Test- 
Model  philosophy  and  practice  promise  to  enlighten  decision-makers’  judgments  concerning  the 
new  system’s  likely  operational  contribution.  This  report  describes  such  an  approach,  and 
quantitative  results  from  them. 
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1. 


Problem  Setting 


The  purpose  of  this  paper  is  to  formulate  a  model  for  the  operational  analysis  of  a 
group  of  mobile  amphibious  vehicles  that  co-operate  on-land  in  a  remote  location,  having 
been  launched  from  the  sea  (littoral/ocean).  The  specific  application  is  to  preview  and 
extend  or  enhance  Operational  Test  and  Evaluation  of  the  Marine  Corps’  Advanced 
Amphibious  Assault  Vehicle  (AAAV),  but  the  model  capability  and  usefulness  extends 
beyond  that  objective.  In  particular,  the  model  can  be  used  to  analyze  logistical  and 
maintenance  support  required  to  carry  out  an  operation  in  a  timely  manner.  The  model  is 
implemented  in  software,  obtainable  from  the  authors,  which  is  executable  on  a  personal 
computer.  The  software  is  used  to  illuminate  logistical  and  maintenance  issues. 

Suppose  a  platoon  of  n  (e.g.,  n  =  12)  nominally  identical  amphibious  vehicles 
(AV  here  an  abbreviation  for  the  AAAV,  but  also  for  other  vehicle  types),  begins  at  time 
/  =  0,  to  traverse  an  ocean  or  littoral  environment  towards  a  beach/shoreline  (B)  a 
distance  d  (miles)  away  from  a  mother  ship  (MS)  source.  Operating  at  high  speed,  the 
vehicles  are  designed  to  travel  at  speed  Vj  (in  miles  per  hour,  which  may  vary,  depending 
on  sea  state,  etc.).  The  ideal  is  that  they  reach  their  first  destination,  B,  without 
interruption;  this  ideal  is  unrealistic,  so  there  is  subsequent  discussion  of  failures  and 
recovery  during  transit  and  their  impact  on  arrival  time  at  the  B.  At  this  point,  they 
prepare  to  become  land  vehicles  by  activating  a  tread  system  and  effectively  becoming 
the  rough,  functional  equivalent  of  Army  Armored  Personnel  Carriers  (APCs).  Their 
individual  payloads  are  platoons  of  Marines.  The  AVs  are  of  smaller  capacity  than  APCs, 
and  are  less  well  armored,  but  carry  more  sophisticated  payloads  (sensors, 
communications,  and  weapons).  They  carry  out  operational  assignments  on  land,  and 
eventually  return  to  a  destination  and  begin  again.  See  Buckles  (1999)  for  further 
discussion  and  references. 
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1.1  Reliability  Considerations 

Realistically,  the  AV  is  made  up  of  failure-prone  subsystems;  propulsion,  drive  train, 
navigation,  communications,  weapons,  etc.;  and  each  subsystem  may  well  contain  one  or 
more  design  faults  that  cause  failure.  Failures  of  any  of  these  subsystems  in  transit 
degrades  the  mission  capability  of  the  AV,  and  hence  of  the  AV  platoon. 

The  impact  of  some  failure  modes  differs  from  others:  certain  failures  of  the  engine 
or  drive  train,  or  other  major  propulsion  system  elements,  incapacitate  a  platform, 
rendering  it  immobile  or  “quiescenf’  (“dead  in  the  water”),  and  require  the  platform’s 
removal  to  a  maintenance  facility.  The  Marine  passengers  are  also  transferred  off  during 
such  a  phase.  Other  failures,  for  instance  of  a  sensor  or  part  of  the  navigation  system, 
may  allow  the  platform  to  keep  up  with  the  remainder  of  the  platoon,  but  in  less  than 
fully  capable  condition.  A  quiescent  period  may  end  after  a  delay  when  an  auxiliary 
vehicle  delivers  spare  parts.  The  present  sophisticated  design  of  an  AV  platform  appears 
to  make  a  great  many  different  failures  potentially  possible,  or  even  likely,  and  to  require 
off-platform  assistance  for  rectifying  such  failures.  Thus,  the  failure-recovery  subsystem 
cannot  be  ignored,  and  should  itself  be  tested  as  part  of  the  overall  system.  Runs  of  the 
proposed  models  will  help  to  indicate  the  effects  of  unreliability  upon  overall  end-to-end 
operational  availability.  Operational  data  acquisition  during  Operational  Testing  (OT)  can 
thus  be  suggested. 

1.2  Purpose  of  the  Paper 

The  major  objective  of  this  paper  is  to  reveal  the  sensitivities  of  overall  AV  system 
operational  effectiveness  to  failure  events  and  alternative  concepts  of  a  recovery-repair 
operation.  It  will  be  seen  that  the  basic  system  of  platforms,  a  platoon,  is  one  of 
potentially  mutually  dependent  subsystems,  the  platforms  themselves  plus  auxiliary 
support.  It  can  also  be  seen  that  if  the  basic  system  is  realistically  expanded  to  include 
support  agents  such  as  landing  craft  and  helicopters,  that  themselves  have  competing 
operational  functions,  then  the  major  surface  and  airborne  subsystems  can  exhibit  strong. 
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interrelated  dynamics  caused  by  reliability-maintainability  system  properties  and 
requirements.  Thus,  careful  acquisition  and  analysis  of  OT  data  is  essential  to  show  that 
there  is  an  acceptable  overall  system  performance  when  such  an  interdependent  system  is 
finally  delivered  to  the  ultimate  users. 

Since  it  is  unlikely  that  an  entire  platoon  of  AVs  and  its  support  system  will  be 
operationally  tested,  it  is  important  that  evaluation  be  supported  by  comprehensive 
and  adequate  system-levels  models,  and  that  these  models  be  parameterized  in  a 
manner  consistent  with  field  experience. 

2.  Model  Outline 

Preliminary  simplified  models  are  presented,  in  which  numbers  of  individual 
platforms  in  various  conditions,  so-called  states,  are  represented  as  functions  of  each 
other,  and  of  elapsed  time  measured  from  platoon  laimch  from  the  MS.  The  basic 
structure  of  the  initial  model  (D-1)  is  chosen  to  illustrate  a  plausible,  dynamic  version 
that  can  be  interpreted  as  a  rough  “expected  value”  (or  deterministic  fluid)  model.  A 
“better”  model  would  be  stochastic,  i.e.,  crafted  to  explicitly  represent  variability  in,  say, 
the  times  between  platform  failures  and  times  to  repair  failed  platforms.  A  simplified, 
analytically  tractable  stochastic  model  (S-1)  is  provided  here.  A  related,  more  detailed 
model,  to  be  investigated  by  computer  simulation  (Monte  Carlo),  will  be  the  subject  of  a 
Modeling  and  Simulation  (M&S)  OR  thesis  by  Capt  Jesse  A.  Kemp,  USMC. 

Next  follows  a  flow  chart  for  the  transitions  of  platforms/vehicles  among  their  various 
states,  as  recognized  in  the  quite  simplified  models  (D-1)  and  (S-1).  Note:  In  (D-1)  we 
ignore  the  effects  of  explicit  distance  or  range  from  the  MS.  Spatial  considerations  will 
be  brought  in  subsequently,  and  are  explicitly  present  in  the  Kemp  thesis. 
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Flow  Chart 


2.2  State  Variables 

A(t):  Number  of  active,  mission-capable  vehicles  in  the  ocean-littoral 

transitioning  at  high  surface  (planing)  speed  across  water  at  time  t. 

L(t):  Number  of  unavailable  (i.e.,  in  transport  service)  auxiliary  vehicles, 

e.g.,  landing  craft  of  type  LCU(R)  (here  abbreviated  LC),  in  the  Littoral 
ocean  at  time  t.  Note:  Available  LCs  number  L  -  L(t)  for  the  present,  L 
being  the  total  number  available.  There  are  other  types  of  LC — the 
LCAs — that  are  not  capable  of  towing  broken-down  AVs.  They  are  not 
considered  at  this  point. 

Qiit):  Number  of  quiescent  AVs  (e.g.,  immobile  or  “dead  in  the  water”),  and 

those  not  fully  mission-capable  at  time  t.  These  are  destined  to  be 
transported,  or  otherwise  proceed  at  low  speed,  to  either  the  MS  or  the  B 
for  repair. 

RaO):  Number  of  AVs  under  transport  by  a  companion  AV  towards  the  MS 

(“in  reverse”)  at  /,  (note  that  RA(t)  counts  the  mxmber  of  vehicles  in  pairs 
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(the  tow-er  “alive,”  the  tow-ee  “dead”)).  Also  note  that  this  number  can 
always  be  zero  if,  by  edict,  A  Vs  do  not  tow. 

Rft):  Number  of  AVs  under  transport  by  an  auxiliary  agent  (e.g.,  landing  craft, 

such  as  an  LCU(R))  towards  the  MS. 

FA(t):  Number  of  AVs  under  transport  (here  tow)  by  another  AV  towards  the  B 

(“forward”)  at  /  (again,  this  counts  individuals  in  pairs).  Note  that  this 
number  can  be  zero.  AVs  may  not  tow,  or  may  do  so  only  if  LCs  (LCUs  or 
LCU(R)s)  are  all  busy  transporting  other  failed  AVs,  or  are  engaging  in 
other  activities. 

FlO):  Number  of  AVs  under  transport  by  an  auxiliary  agent  (e.g.,  landing  craft) 

towards  the  B. 

Ms(t):  Number  of  AVs  undergoing  or  awaiting  (queuing  for)  maintenance  on  the 

MS  at  t. 

Msd):  Number  of  AVs  undergoing  or  awaiting  maintenance  on  the  B  at  t. 

B(t):  Number  of  operational  AV  platforms  accumulated  on  the  B  at  t.  Some 

have  arrived  without  failure  (perhaps  after  being  towed  to  the  MS, 
repaired  and  relaunched);  some  have  experienced  maintenance  on  the  B 
after  arriving  there. 

2.3  Parameters 

Here  are  parameters  minimally  describing  the  dynamic  evolution  of  the  system  of 
platforms,  as  it  transitions  (planes)  from  the  MS  to  the  B. 

A,.-  Failure  rate  per  platform  while  running  at  fast  speed;  1/  A,  is  the  “mean 

time  to  failure”  (MTTF),  measured  from  any  time  point  at  which  a 
platform  is  up  and  running.  (Units  are  1/time  =  e.g.,  1/hours).  Note:  The 
value  of  X  may  depend  upon  speed,  sea  state,  condition  of  the  littoral 
ocean  environment,  and  other  explanatory  variables;  this  can  be  written  as 
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A(»).  Further,  a  model  that  permits  infant  mortality  (early  post-launch 

failure)  can  easily  be  adapted  from  the  present  setup.  Exercise  of  the 
resulting  model  will  show  the  possible  considerable  sensitivity  of 
outcomes  to  a  natural  departure  from  the  usual  constant-failure-rate 
(exponential)  assumption. 

Pi:  Probability  of  “infant  failure,”  meaning  that  a  platform  fails  after  launch, 

to  attain  planing  speed.  It,  therefore,  immediately  joins  a  queue  of  AVs 
awaiting  maintenance  at  the  MS. 

vpa:  Rate  at  which  quiescent/failed  AVs  are  picked  up  (placed  under  transport, 

i.e.,  towed)  by  fellow  AVs  (this  can  be  omitted/made  zero,  thus 
eliminating  the  AV-tow  option  by  setting  vpa  =  0).  (Units  are  1/time  = 
e.g.,  1 /hours). 

vpp:  Rate  at  which  quiescent/failed  vehicles  are  picked  up  (placed  under 

transport,  i.e.,  towed)  by  auxiliary  vehicles  (e.g.,  LCs,  such  as 
the  LCU(R)). 

v/>  =  vpA  +  vpi:  Total  rate  of  tow-er  pickup  by  tow-ee. 

Vf.'  Rate  at  which  individual  platforms  reach  the  B;  1/vr  is  the  mean 

time  to  transition  to  the  B  while  in  planing  motion  (given  that  no 
failure  occurs).  (Units  are  1/time  =  e.g.,  1 /hours),  vj  can  be  set 
equal  to  V-fd  to  calibrate  the  present  model  to  one  with 
deterministic  motion,  distance  to  the  beach  d,  speed  velocity  Vp, 
and  no  failures  (see  Section  4  for  related  discussion). 

Ps(*):  Rate  at  which  the  MS  maintenance  facilities  can  return  a  failed 

platform  to  fully  active  capability. 

Note  1:  For  the  present,  MS  maintenance  is  treated  as  a  saturable, 
single-server  repair  system. 
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Par(*): 

Pab(*): 

Plr(*): 

Plb('): 


Note  2:  When  there  is  infant  mortality,  i.e.,  pi>0,  and  a  platform 

fails  essentially  just  after  being  launched,  then  the  effective  repair 
rate,  equals  an  actual  repair  rate  multiplied  by  the 

probability  that  no  infant  failure  occurs: 

Ps[*)  =  {^- Pi)p%*)-, 

in  other  words,  the  effective  mean  repair  time  (effective  repair  time 
random  variable  is  denoted  by  Rgff)  with  infant  failures  present, 

with  probability  is 

This  is  because  an  average  of  l/(l-j9,)  repair  times  must  be 

performed  before  the  infant-failure  mode  is  avoided  (the  present 
model  does  not  allow  for  that  mode’s  permanent  removal;  see 
Gaver  et  al.  (2000)  for  plausible  reliability  growth  alternatives). 
Further  Note:  The  notation  refers  to  the  influence  of  other 
variables,  such  as  feedback,  to  enhance  service  if  a  big  backlog 
develops.  The  default  in  all  cases  is  a  constant  parameter  value. 
Maintenance  rate  at  the  B. 

Probability  that  a  failed  vehicle  is  towed  to  the  MS  (reverse)  by  a 
fellow  AV  (two  AVs  in  complex,  one  “alive”  and  one  “dead”). 
Probability  that  a  failed  vehicle  is  towed  to  the  B  by  a  fellow  AV. 
Note:  par(*)  +  Pab(')=  1- 

Probability  that  a  failed  vehicle  is  transported  to  the  MS  (reverse) 
by  an  auxiliary  vehicle  (e.g.,  LC). 

Probability  that  a  failed  vehicle  is  transported  to  the  B  by  an 
auxiliary  vehicle  (e.g.,  LC).  Note:  Plr(*)  +  Plb(*)=  1- 
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v^;  Rate  at  which  a  towed  vehicle  complex  (two  A  Vs,  one  “alive”  and 

one  “dead”)  reaches  the  MS;  1/  v/m  is  the  mean  time  for  the  towed 
complex  to  reach  the  MS. 

vri:  Same,  when  transport  is  by  auxiliary  vehicles  (e.g.,  LC). 

vba-'  Rate  at  which  a  towed  vehicle  complex  reaches  the  B;  1/  vba  is  the 

mean  time  for  the  towed  complex  to  reach  the  B. 
vbl-'  Same,  when  transport  is  by  auxiliary  vehicles  (e.g.,  landing  craft). 

L:  Total  number  of  landing  craft  used  for  AV  recovery  (LCU(R)s,  not 

to  be  confused  with  LCACs)  available  in  the  ocean  littoral  (assume 
fixed  for  present). 

Simplification  1.  Pyramiding/cascading  failure  possibilities  are  not  modeled  here.  If 
a  failed  and  quiescent  vehicle  must  be  transported  or  towed  without  failure,  it  is 
presumed  transported  or  towed  to  the  MS  or  to  the  B,  depending  upon  which  destination 
is  perceived  as  the  best  from  the  viewpoint  of  total  vehicle  soonest  availability  at  the  B 
waypoint.  The  problem  is  potentially  made  more  complex  by  the  possibility  that 
previously  failed  towed  vehicle  complexes  (e.g.,  two  AVs,  or  an  AV  and  LC,  may  be 
ahead  of  a  new  breakdown,  and  hence  virtually  queued  up  ahead  of  that  breakdown  in 
time).  Consequently,  a  breakdown  that  occurs  near  the  MS  might  still  be  advantageously 
towed  to  the  B  if  there  are  other  (higher  priority)  transported  or  towed  complexes  between 
the  current  breakdown  and  the  MS. 

Simplification  2.  It  is  assumed  here  that  a  transported  or  towed  complex  does  not 
break  down  between  its  origin  point  and  destination  (MS  or  B).  This  is  optimistic,  and 
calls  for  future  amendment,  although  the  model  then  becomes  more  complex  and  difficult 
to  manage.  Failures  that  are  not  total  breakdowns  can  be  handled  in  terms  of  the 
current  model. 


10 


Simplification  3.  For  the  purposes  of  this  discussion  it  will  simply  be  assumed  that  if 
a  breakdown/failure  occurs  it  is  towed  to  the  MS  with  probability  p.R,  and  to  the  B  with 
probability  p.B,  where  refers  to  the  tow-er  type  Par(*),  Plb(*)  (see  2.3).  A 
convenient,  but  very  preliminary  assumption  is  that  these  probabilities  are  1/2.  However, 
it  is  clear  that  more  realistic  probabilities  should  depend  minimally  on  the  time  that  has 
elapsed  since  platoon  launch  from  the  MS;  the  probabilities  can  also  be  made  to  depend 
upon  repair  backlogs  and  delays,  possibly  due  to  past  availability  at  the  two  possible 
destinations.  In  general,  a  more  refined  description  of  the  system  state  should  lead  to  a 
better  decision,  but  at  a  cost  in  model  complexity. 

Clearly  the  decision  as  to  where  to  tow  the  failed  AV  should  ultimately  depend  upon 
how  soon  it  can  be  in  service  near  (e.g.,  inland  of)  the  B. 

Simplification  4.  The  mean  time  to  be  towed  to  either  the  MS  or  the  B  should  depend 
systematically  and  predictably  upon  the  actual  breakdown  location,  although  in  practice 
there  may  be  great  variability.  As  a  surrogate,  the  mean  time  could  depend  on  elapsed 
time,  t;  presumably  as  elapsed  time  t  increases,  the  distance  of  the  platoon  from  B 
decreases  {but  not  if  early  breakdowns  persistently  occur  near  the  MS,  since  the  repair 
facility  there  may  become  overloaded).  Our  present  assumption  of  a  constant  mean  time 
is  thus  crude,  since  it  merely  represents  an  order-of-magnitude  delay.  The  assumption  can 
be  relaxed  by  segmenting  the  Littoral  region  into  several  range  bands  from  the  MS  to  the 
B  (and  beyond),  and  counting  platforms  in  various  states  in  each.  This  refinement 
proliferates  the  state  equations,  and  is  left  for  a  second  round  of  modeling.  It  will 
subsequently  be  shovra  that  the  simple  Markovian  type  of  fluid  model  presented  here,  is 
in  quite  good  agreement  with  a  more  limited,  but  stochastic,  model. 

Many  of  the  simplifications  mentioned  have  been  avoided  in  a  Monte  Carlo  computer 
simulation  model  implemented  by  Capt  Jesse.  A.  Kemp  (USMC)  as  part  of  his  Masters 
thesis  in  Operations  Research. 
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3.  Deterministic  Analytical  (Mathematical)  Model  (D-1) 

Here  is  a  set  of  dynamic  stock-and-flow  differential  equations  that  describe  the 
time-dependent  numbers  of  AV  platforms  in  the  various  states.  The  present  model  “stops 
at  the  beach,”  i.e.,  at  the  B.  The  model  is  continued/amended  in  Section  4,  to  include  the 
on-land  mission  segment. 


^  =-^(0 

V - - / 

'' — V — ^  Rate  of  failure 

Rate  of  "increase"  of  active  AVs 

of  number  of  active 
AVs  in  the  Littoral 


''  V 

Rate  of  tow  initiation 
for  quiescent  (broken  - 
down)  AVs 


Rate  of  "free"  active 
AVs  reaching  B 


(3.1) 


-h 


Ms{t) 


Rate  at  which  AVs  under 
maintenance  at  MS  become 
active  in  the  Littoral 


+  ^ra^aA) 

— 

Rate  at  which  AVs 
under  tow  by  fellow 
AVs  reach  the  MS 
(reverse  path);  tow-er 
released  immediately 


dQT(t) 

Rate  of  "increase"  of 
number  of  quiescent 
AVs  awaiting  tow 


Am  -  (t) A(t)  -  (t)[L  -  L(t)f 

Rate  of  towable-  Rate  of  tow  initiation  Rate  of  removal  initiation 

failure  events  by  active  AVs  by  available  auxiliary 

vehicles  (e.g.,  LCs) 


(3.2) 


Note:  L(t)  =  RlO)  +  Fi(t).  Also,  0  <L(t)  <L, 

so  [L-L(t)f=  maxfL  -  L(t),  0] 

Note:  Suitable  initial  conditions  are 
A{Q)  =  {\-p)n 

Mv(0)  =  p^n 


where  n  is  the  number  of  AVs  launched  from  the  MS 


—4— 

Rate  of  "increase"  of 
number  of  AVs 
(tow-er  and  tow-ce) 
being  towed  ("reverse") 
to  the  MS 


2PAR(*)VpAQT(t)A(t) 

Rate  of  initiation  of 

of  tow  (complex  of  two  vehicles) 

towards  the  MS 


VraRaW 

\ - - / 

Rate  of  dropoff  of 
AVs  for  maintenance 
at  the  MS 


VraRaW 

Rate  of  tow-er 
return  to  active 
status,  and  to  hover 
rate  towards  the  B 


(3-3) 


Note:  When  an  active  vehicle  picks  up  and  tows  a  failed,  quiescent  vehicle,  this  means 
that  two  vehicles  (AVs)  form  a  complex  and  proceed  at  rate  vpa.  Hence,  the  “2.” 
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=PLK(»)Vp,.QT(t)[L-L(t)]- 

^  '  Rate  of  pickup  of  failed 

Rate  of  ”  mcrease”  of  A  Vs  by  available  landing 

the  number  of  AVs  craft:  headed  for  the  MS 

being  transported  to 
the  MS  by  an  auxiliary 
vehicle 


VrlRlO) 

V _ _ * 

Rate  of  delivery 
of  landing  craft 
towed  AVs  to 
the  MS 


(3.4) 


Ml 

Rate  of  "increase" 
of  number  of  AVs 
in  tow  ("forward") 
to  the  B 


=  2p^(»)vpAQT(t)A(t)  -  VBApACt) 

\ ^ ^  V  ^ 

Rate  of  initiation  of  Rate  of  dropoff  for 

tow  (complex  of  2  vehicles)  maintenance  at  the 

towards  the  B  B 


'■  V . . . 

Rate  of  tow-er 
return  to  active 
status 


(3.5) 


'Note:  In  this  case,  both  tow-er  and  tow-ee  remain  at  the  B,  but  the  tow-ee  may  enter 
repair,  encountering  delay.  counts  both  tow-ers  and  tow-ees. 


=  PLni^ypLQA^tL-  L{t)]\ 

'' — '  Rate  of  initiation  of  transport  by  LC 

Rate  of  "increase"  of 
number  of  AVs  in  tow 
to  B  by  LC 


- - - / 

Rate  of  dropoff 
of  LC  tows  at  B 


(3.6) 


“'^RA^aCO  +  ^rl^l(^) 


As(«)Ms(t) 

1  +  Ms(t) 

' - - - ' 


dM,(l) 

dt 

Rate  of  "increase" 
of  number  of  AVs 
in  MS  maintenance 


- v - 

Rate  of  dropoff 
for  maintenance 
at  the  MS  from 
AV  tow 


V 

Rate  of  dropoff  at 
MS  from  LC 
transport 


Rate  of  maintenance 
completion  at  the  MS 


(3.7) 


dMB(t) 

Rate  of  "increase" 
of  number  of  AVs 
in  B  maintenance 


^ba^aW 

V - ^ / 

Rate  of  AV  dropoff 
for  maintenance 
at  the  B  by  AV 
tow-ers 


^BLFL(t) 

' - ^ ' 

Rate  of  AV  dropoff 
for  maintenance  at 
the  B  by  LC 
tow-ers 


/^B(*)MB(t) 

l  +  MB(t) 

_  -  ^ 

V" 

Rate  of  maintenance 
completion  at  the  B 


(3.8) 


-4- 

Rate  of  increase 
of  LCs  unavailable 
for  failed  AV 
transport  (in  use) 


Rate  of  dropoff 
of  AVs  by  LCs 
at  the  MS 


Rate  of  dropoff  R^te  of  pickup  of  quiescent /failed  AVs 
of  AVs  by  LCs  by  available  LCs 
at  the  B 


(3.9) 
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Rate  of  "increase" 
of  active,  available 
AVs  at  the  B 


=  MW  +,V^BAFA(t)  + 

Rate  of  "free"  Rate  of  active 

active  AVs  tow-er  AVs 

reaching  the  reaching  the 

B  B 


I  +  MbW 

- 

Rate  of  maintenance 
completion  of  the  B 


(3.10) 


4.  The  On-Land  Mission.  Helicopters  (For  Transport  and  AV 
Rescue-Recovery).  Model  D-2. 

An  important  aspect  of  Marine-Navy  cooperation  to  dominate  and  secure  a  land 
region  “from  the  sea”  is  (in  present  concept)  to  use  AVs  and  Helicopters  (Hs  for  short), 
or  some  other  mobile  platform,  cooperatively.  A  primary  function  of  an  H  force  is  to 
transport  personnel  and  supplies  to  inland  points.  An  important  secondary  role  is  to 
support  AVs,  when  the  latter  fail,  while  in  the  on-land  mission  phase. 

We  now  aim  to  more  explicitly  include  in  the  model  the  contribution  of  Hs  to  AV 
dynamics,  and  ultimately  to  AV  availability  on  the  B.  Conversely,  we  also  study  the 
impact  of  AV  requirements  for  support  on  H  availability.  There  will  clearly  be  a  tradeoff 
that  can  be  dependent  upon  conflict  conditions  ashore.  The  present  model  does  not  reflect 
the  time  and  event-dependent  character  of  such  needs,  but  only  general  trends  that  result 
from  a  consistent  control  policy. 

4.1  Augmentation  of  the  Basic  Model  to  Represent  AV  Operational  Region 
(On-Land)  Availability 

The  focus  here  is  on  AVs  that  may  fail  when  on  land,  and  can  be  repaired  on-site. 
This  accounts,  at  least  crudely,  for  the  use  of  Hs  to  furnish  needed  parts  to  a  stranded  AV. 

Now  consider  the  availability  of  AVs  that  have  eventually  and  successfully  reached 
the  B,  and  undergone  successful  repair,  if  needed.  The  state  variable  B(t)  of  (3.10)  counts 
that  number  at  time  t,  provided  none  are  yet  dispatched  to  carry  actual  in-land  operational 
missions.  But,  such  are  the  ultimate  objectives  of  the  AVs,  so  we  now  add  features  to  the 
model  to  represent  this  phase  of  the  AV  operational  cycle. 
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4.2  Additional  State  Variables 

These  augment  those  of  Section  2.2  to  represent  on-land  conditions  (the  subscript  O 
denotes  this  on-land  operational  phase  throughout); 

Ao(t):  Number  of  active,  mission-capable  vehicles  (AVs)  on  land,  transitioning 

among  waypoints  at  time  t.  These  are  assumed  to  be  failure-prone  with 
rate  ■ 

Oo(t):  Number  of  failed,  quiescent  AVs  (e.g.,  in  place  on  land)  and  not  fiilly 

mission  capable  at  t.  These  require  recovery  assistance;  we  assume  it  is 
furnished  on  site  by  helicopters  or  other  support  vehicles.  Other  vehicles 
(e.g.,  for  transport)  can  also  be  accommodated  in  the  model,  but  are 
not  here. 

Let  there  be  h  (e.g.,  five  or  10)  helicopters  (Hs)  or  other  mobile  support  vehicles 
exclusively  dedicated  to  assist  failed  on-land  AVs. 

Let 

Ho(t):  Number  of  Hs  outgoing  to  assist  failed  AVs  at  t. 

Hi(t):  Number  of  Hs  incoming  from  having  assisted  an  on-the-spot  repairable 

AV,  at  t. 

Note:  It  is  assumed  here  that  H  transit  times,  both  outbound  and  incoming,  are  much 
greater  than  the  time  on  site,  which  is  taken  to  be  negligible. 

h-Ho(t)-Hi(t):  Number  of  Hs  available  to  be  sent  on  AV  recovery/assistance  missions, 
at  t. 
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4.3  Parameters 

Let 

/Iq(*);  The  failure  rate  of  AVs  on  land.  As  before,  the  notation 

stands  for  the  influence  of  explanatory  variables,  such  as  terrain 

state  and  speed  of  advance.  Here  the  default  is  a  constant. 

(•).•  Rate  at  which  outbound  Hs  encounter  failed/quiescent  AVs 

on  land. 

Vh  (•).'  Rate  at  which  inbound  Hs  return  to  H  pool,  and  become  available 

for  reassignment  to  AV  recovery. 

Rate  at  which  active  AVs  leave  mission  region  (return  to  the  B  or 

reach  a  destination);  1  /  v/,  is  the  mean  mission  duration. 

Note:  The  h  Hs  are  considered  to  be  exclusively  assigned  to  provide  on-land,  on-site 
recovery  assistance  to  failed  AVs.  They  are  themselves  assumed  immune  to  crashes  and 
failures  (an  unrealistically  optimistic  assumption  in  hostile  territory,  but  one  that  can  be 
relaxed).  They  are  also  assumed  to  make  only  one  trip — ^and  that  successful — ^per 
failed  AV. 

Let 

0  <  r'  <  tj:  (Denotes)  the  rate  at  which  available  AVs  on  the  B  enter  on-land 
mission  status  at  time  t. 

Here  are  state  equations  to  describe  the  on-land  mission  availability  of  the  AVs. 


=  f?7(B(t'),  0<t'<t)  -  Ao(*)Ao(t)  -  viAo(t) 

' - '  Rate  of  mission  force  creation  Rate  of  failure  of  Rate  of  depletion 

from  those  on  B  active  on~land  AVs  of  AVs  in 

mission  region 


(4.1) 


Rate  of  H -assisted  recovery  of 
failed-quiescent  AVs 
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(4.2) 


^=-  y«A’‘)Qo{iW') 

at  ^ - V - '  ' - V - ' 

'  '  Rate  of  H-assisted  recovery  of  Rate  of  failure  of 

failed -quiescent  AVs  active  AVs 


dHo{t) 

dt 


VHo  {•)Qo{t)Ho{t)  +  Ao{*)Ao{t)l[h  -  Ho{t)  - 

' - V— ;; - '  ' - V - ^ 

Rate  of  H-assisted  Rate  of  launch  of  Hs  in  response  to  failures 

recovery  of  failed  AVs 
(rate  of  outgoing  H 
release) 


(4.3) 


Note: 


That  is,  an  H  recovery  mission  begins  only  if  at  least  one  or  more  Hs  are  available. 


(4.3a) 


^  ^  Rate  of  H-assisied 

recovery  of  failed  AVs 
(rate  of  outgoing  H 
release) 


V - ^ - / 

Rate  of  return  of  Hs 
to  base 


(4.4) 


4.4  AV  Beach  Accumulation  and  Land  Mission 

Although  the  platoon  of  n  (e.g.,  n=12)  AVs  departs  from  the  MS  essentially  as  a 
concentrated  group,  its  members  accumulate  gradually  at  the  B.  In  the  fluid 
approximation  model  used  here,  the  accumulated  number  on  the  B  approaches,  but  never 
actually  reaches,  «;  this  is  because  a  continuous  “fluid”  approximation  is  being  used.  We 
adopt  the  convention  that  “all”  AVs  are  effectively  at  the  B  when  B(t)  reaches  nfm,  where 
0  <  ^  1 ,  but  /„,  is  taken  to  be  large,  e.g.,  0.95  or  0.99  (choice  of  the  analyst).  At  this 

point,  the  on-land  platoon-strength  mission  is  launched. 

Thus,  (a)  calculate  the  launch  time,  : 

^{fm  )  ~  hfm  » 

then  (b)  begin  integrating  (4.1)  at  r  =  and  invoke  the  initial  condition  Ao{tm)  =  so 
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Initial  Conditions  (I) 


0<t<t^ 

t>tm 

o 

11 

Ao{t)  =  nf„ 

Solve  differential  equations 

o 

11 

Qo[t)  =  0 

o 

11 

Ho{t)  =  0 

4.5  Upper  Bound  on  Number  Available  During  On-Land  Mission  Segment 

In  order  to  obtain  a  deterministic  upper  bound  on  the  on-land  availability  of  the 
platoon,  pick  /„,  =  1  - 1  /  2« ,  which  stops  B(t)  accumulation  when  B{t^^  —  n-\l2,  an 

approximation  to  the  mean  first-passage  time  to  n.  Thus,  (a)  launch  a  mission  at 
where 

=  =«-l/2, 

and  use  the  initial  condition  AQ{t^)  =  n.  Then,  (b)  solve  (4.1)  -  (4.4).  These  define 
Initial  Conditions  (II). 

The  two  simple  initial  conditions  specified  in  Subsections  4.4  and  4.5,  allow  the 
general  term  in  (5(C),  0  <  t'  <  /)  in  (4.1)  to  be  omitted;  integration  starts  at  time  with  a 

given  level  of  platoon  strength.  Note:  The  above  specifies  a  particular  and  plausible 
initial  condition  for  the  on-hand  mission  segment  of  the  AV  platoon.  Other  modes  of 
dispatch  of  the  A  Vs  into  the  on-land  environment  can  be  represented  by 
straightforward  changes. 
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4.6  Numerical  Illustration  of  the  Sensitivity  of  Operational  Availability  to  the 
Form/Shape  of  the  Time-to-Failure  Distribution 

In  this  subsection  we  compare  results  for  an  exponential  time  to  failure  to  those  of  a 
mixed  failure  time  distribution,  in  which  there  is  a  probability  /?,  that  the  failure  time  is  0 
and  a  probability  (1-p,)  that  the  failure  time  is  exponential.  The  exponential  and  the 
mixed  exponential  have  the  same  means,  but  their  shapes  are  different.  The  deterministic 
model  of  Sections  3  and  4  is  used.  There  are  12  AVs  in  a  platoon.  Model  parameters 
common  to  all  the  cases  are  as  follows: 

AVs  traveling  to  the  B: 

Distance  from  the  MS  to  the  B  =  40  nautical  miles  (nm) 

All  failed  AVs  in  water  are  towed  back  to  the  MS 

AV  velocity  in  water  =  25nm/hr 

Towing  velocity  in  water  =  Snm/hr 

If  an  AV  does  not  fail,  the  time  to  reach  the  B  is  1 .6  hours 

Arrival  rate  at  the  B,  vy-  =  25  /  40 = 0.625 

Towing  rate  back  to  the  MS  if  failed,  v/>^  =  5/20  =  0.25 

Stop  at  the  B  until  10  of  the  AVs  are  available  at  the  B;  then  the  available  AVs  start 
the  land  portion  of  mission 
AVs  on  Land: 

Distance  from  the  B  to  last  destination  =  1  OOnm 

AV  velocity  on  land  =  25nm/hr 

Travel  time  on  land  without  failure  =  4  hours 

There  are  five  waypoints;  an  AV  pauses  10  minutes  at  each  waypoint 

Total  time  spent  at  waypoints  =  50  min. 

AVs  are  subject  to  failure  at  the  waypoints 

Total  time  an  AV  transits  on  land  if  no  failure  =  4.83  hours 

Rate  at  which  an  AV  arrives  at  land  destination,  vy  =  1/4.83  =  0.21 
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Helicopter: 

Velocity  =  SOnm/hour 

Mean  time  to  fly  to  failed  AV  =  (40+(100/2))/50  =  1.8  hours 
There  is  a  two-hour  administrative  logistic  down  time  (ALDT)  prior  to  each  H 
mission  to  repair  a  failed  AV  on  land 

Rate  of  arrival  of  an  H  from  the  MS  to  a  failed  AV,  vh.  =  l/(2+l  .8)  =  0.26 
Rate  of  return  of  an  H  to  the  MS,  vh„-  1/1-8  =  0.56 
There  are  two  Hs 

Results  for  the  following  AV  time  to  failure  models  are  compared;  all  models  for  time 
to  failure  for  an  AV  have  a  mean  of  72  hours.  Two  of  the  failure-time  models  have  infant 
failure  modes  as  described  in  Section  3. 

Exponential  model: 

Mean  time  to  failure=72  hours 
Mean  repair  time  on  the  MS=three  hours 
Failure  rate,  A=l/72=0.014 
Rate  of  repair  completion  on  the  MS=l/3=0.33 
Initial  condition:  all  12  AVs  are  available  at  time  0 
Mixed  model  1: 

With  Pi  =  0.25  the  failure  time  is  0,  and  with  probability  (7-/?/)  =  0.75  the  time  to 
failure  is  exponential  with  mean  96  hours;  the  mean  time  to  failure  is  72  hours 
Failure  rate,  X=\I96  =  0.010 

Rate  of  repair  completion  on  MS  (includes  repair  extensions  due  to  additional  failures 
at  time  0)  =  0.75/3  =  0.25 

Initial  condition:  nine  AVs  are  available  at  time  0,  and  three  are  waiting  for  or  are 
being  repaired. 
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Mixed  model  2: 


With  Pi  =  1/12  the  failure  time  is  0,  and  with  probability  (/-/?,)  =  1 1/12  the  time  to 
failure  is  exponential  with  mean  78.5;  the  mean  time  to  failure  is  72  hours. 

Rate  of  repair  completion  on  MS  (includes  repair  extensions  due  to  additional  failures 
at  time  0)  =  (ll/12)/3  =  0.31 

Initial  condition:  1 1  AVs  are  available  at  time  0,  and  one  is  being  repaired. 

The  following  figures  display  the  accumulated  mean  number  of  AVs  that  arrive  at  the 
B;  the  mean  number  of  AVs  in  repair  at  the  MS;  and  the  mean  number  of  AVs  in  transit 
on  land  to  the  final  land  destination. 

Note:  The  form  of  the  failure  distribution  greatly  influences  the  time  at  which  the 
platoon  reaches  the  final  land  destination.  The  assumption  of  a  “pure  exponential” 
distribution  tends  to  be  quite  optimistic  if  an  “infant  failure”  type  of  alternative  holds. 
Conclusion:  OT  data  acquisition  and  analysis  should  allow  a  propensity  for  early 
failures,  as  compared  to  the  exponential,  to  reveal  itself.  Automatic  assumption  of  a 
simple  exponential  may  be  dangerously  optimistic.  See  Section  6  for  a 
similar  message. 
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Number  of  AVs  to  Arrive  on  the  Beach  (B)  Number  of  A  Vs 


0  20  40  60  80  100  120 

Time  (Hours  from  MS  Launch) 


—  Mixed  dist  with  mean  72;  with  prob  pi=0.25  time  to  fail  -0,  o/w  expon  with  mean  96  j 

—  Exponential  time  to  failure  mean  72  | 

’  Mixed  dist  with  mean  72:  with  prob  pi=l/l  2  time  to  fail  ==0,  o/w  expon  with  mean  78.5  [ 

Figure  1. 

Number  of  AVs  in  Repair  at  the  Mother  Ship  (MS). 


Exponential  time  to  failure  mean  72 

Mixed  dist  with  mean  72:  with  prob  pi=0.25  time  to  fail  -0,  o/w  expon  with  mean  96 
Mixed  dist  with  mean  72:  with  prob  pi=l/I2  time  to  fail  =0,  o/w  expon  with  mean  78.5 


Figure  2. 

Cumulative  Number  of  AVs  that  Arrive  on  the  Beach  (B). 
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i - Mixed  dist  with  mean  72;  with  prob  pi=0.25  time  to  fail  =0,  o/w  expon  with  mean  96  j 

‘ - Exponential  time  to  failure  mean  72  ! 

. . Mixed  dist  with  mean  72:  with  prob  pi=l/12  time  to  fail  -0,  o/w  expon  with  mean  78.5  ; 

Figure  3. 

The  Mean  Number  of  AVs  in  Transit  on  Land. 

5.  A  Stochastic  Model  for  Platform  Transit  Time 

It  is  easily  possible  to  extend  the  previous  analytical  model  to  a  tractable  stochastic 
version  that  allows  for  general  (not  just  Markov/exponential  constant  failure  rate) 
random  times  to  failure.  Restriction  to  the  situation  in  which  individual  platforms  operate 
independently  is  immediately  tractable,  and  quite  informative;  it  shows  quantitatively 
that  the  ultimate  transit  time  to  B,  Td  may  be  quite  sensitive  to  the  form  of  the  time-to- 
failure  distribution.  Facile,  knee-jerk-conventional  assumption  of  “the  exponential”  to 
describe  time  to  breakdown,  Td,  can  considerably  understate  the  (expected)  time  for  a 
platform  to  ultimately  reach  the  B  from  the  MS. 

Consequently,  it  is  important  to  tabulate,  study,  and  incorporate  actual, 
observed  times  to  platform  failure,  particularly  those  that  require  unanticipated 
early  return  to  the  MS,  and  subsequent  repair  and  re-launch.  An  important 
objective  of  testing  is  to  identify  and  remove  the  causes  of  such  so-called  “infant 
failures.” 

It  can  be  seen  intuitively  that  a  particular  mean  time  to  failure  can  be  consistent  with 
(a)  a  great  preponderance  of  short  time  values,  balanced  by  some  (enough — even  one!) 
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veiy  long  observations  (sometimes  called  the  infant  mortality  situation)  and  (b)  an 
apparently  exponential-like  scatter,  or  (c)  numerous  alternatives.  Clearly  option  (a),  with 
its  many  repairs  and  re-starts,  can  radically  extend  the  time  to  waypoint  B  by  inducing 
many  failures  that  must  await  assistance  or  return  to  the  MS,  experience  maintenance  and 
re-launch.  An  important  exception  occurs  if,  after  each  such  “infant”  short-time  failure,  a 
serious  attempt  is  made  to  remove  a  design  fault  that  causes  the  failure,  and  that  attempt 
is  successful  (and  does  not  somehow  otherwise  disrupt  the  design).  This  highly 
reasonable  concept  of  design  evolution  and  testing  (“test-analyze-fix-test”)  can,  it  is 
hoped,  result  in  non-homogeneous  time-to-failure  data  that  exhibits  a  general  increase  in 
operating  times  to  failure  (so-called  “reliability  growth”).  Certainly  the  summary  of  such 
data  by  a  simple  mean,  and  the  invocation  of  a  constant  failure  rate  1  =  1/  mean,  are 

overly  simplistic  and  quite  possibly  misleading.  {Note:  This  indicates  that  our  first  model 
(D-1)  of  Section  3  must  be  viewed  with  caution;  however,  it  can  be  readily  altered  to 
reflect  infant  mortality,  or  relatively  many  premature-compared-to-exponential  failures). 

As  was  the  case  in  Model  (D-1),  some  failures  are  serious,  and  require  slow  transit, 
e.g.,  back  to  the  MS,  or  forward  to  the  B,  for  repair.  This  can  be  accomplished  by  towing 
by  a  fellow  non-disabled  AV,  or  transported  by  LCU(R),  as  represented  in  Model  (D-1) 
in  Section  3.  In  the  present  “independent”  model,  the  failed  AV  is  assumed  to  proceed 
under  its  own  or  totally  independent  power,  or  with  LCU(R)  transport  assistance,  so  the 
AV  platoon  is  not  depleted  for  towing.  This  represents  the  situation  in  which  there  are 
always  ample  LCU(R)s  or  other  vehicles  for  transport  of  failed  AVs.  Its  failure- induced 
slow-retum  velocity  is  Wr  «  Vr.  Minor  (not  serious)  failures  can,  in  principle,  occur 
while  a  seriously  failed  AV  is  being  returned  to  the  MS,  and  thus  adding  the  maintenance 
burden,  and  hence  delay  at  the  MS.  An  interesting  Concept  of  Employment  (COE)  policy 
question  is  whether  to  (a)  service  on-site  by  LCU(R),  minor  on-site  failures  soon  after 
occurrence,  without  transporting  the  AV  to  the  MS  or  the  B,  or  (b)  to  wait  until  the  AV 
reaches  the  MS  or  the  B,  and  service  any  serious  failure  and  the  accumulated  minor 
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casualties  at  that  time.  Of  course,  the  same  issue  arises  if  serious  failures  are  towed  to  the 
B  for  repair:  should  auxiliary,  e.g.,  LCU(R),  support  be  used  to  service  minor  failures  as 
they  occur  (or  en  route),  or  should  the  service  of  accumulated  minor  plus  serious  failures 
be  delayed  until  the  B  is  reached?  This  latter  COE  reduces  the  burden  on  auxiliary 
vehicles,  but  might  jeopardize  partially  disabled  AV  safety,  and  also  potentially  submits 
the  auxiliary  vehicles  to  failure.  The  control  policy  should  depend,  in  part,  upon  current 
need  for  auxiliary  vehicles  for  other  missions,  and  on  the  maintenance  backlog,  hence 
delay  at  the  MS  or  the  B.  Logistics  (spare  parts  and  diagnostics  equipment)  requirements 
must  be  considered  as  well,  and  may  be  addressed  in  terms  of  the  current  models  and 
natural  variations  on  them.  Many  such  are  possible  and  remain  to  be  analyzed. 

5.1  Stochastic  Model  (S-1) 

Here  is  a  simple  stochastic  model  for  recovery  of  a  remotely  failed  AV  during  ocean 
transit.  A  single  generic  platform/vehicle  travels  at  constant  velocity  Vj  (miles/hour)  from 
the  MS  to  the  B,  a  distance  d.  This  distance  can  be  covered  in  precise  time  t  =  d/Vr, 
unless  there  are  interrupts  caused  by  system  failures,  which  we  must  account  for. 
Note:  Here  we  assume  that  actual  time  of  transit  when/if  no  failures  occur  is  a  constant, 
whereas  in  the  previous  model  that  time  is  mean  (of  an  exponential  random  variable). 

Notation 

Td:  Random  time  for  a  representative  platform  to  reach  the  B  from  the  MS  for 

the  first  time,  having  started  at  launch  (time  =  0),  and  possibly  experienced 
various  numbers  of  failures  and  repairs  en  route.  Discovery  of  the 
probabilistic  properties  of  this  random  variable  is  the  objective  of  the 
model.  Model  nms  can  then  guide  analyses  of  alternative  COEs. 

Tb:  The  random  time,  measured  from  a  laimch  from  the  MS  until 

breakdown/failure,  while  attempting  to  transit  to  the  B; 
>t}  =  \- Fj.^{t  =  d  IVj),  where  iv^(*)  is  the  distribution  of  the 
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(arbitrary)  time  to  breakdown  (under  stable  design  conditions),  more 
generally  where  introducing  a  generic  parameter  9 can  represent 

any  quantified  set  of  conditioning  or  explanatory  variables,  or  it  can 
represent  a  random  environmental  effect  as  well.  A  differently  distributed 
Tb  will  prevail  over  land.  The  present  model  pertains  to  the 
ocean-crossing  segment  of  the  transit  only. 

Rq:  Random  time  to  rectify  a  failure  that  is  susceptible  to  on-the-spot  repair, 

given  spare  parts  and  diagnostic  capability.  This  incorporates  the  time  for, 
say,  LC  transit  to  the  failed  AV,  and  even  possibly  necessary  revisits. 

Rm\  Random  time  to  repair  a  serious  failure  that  cannot  be  rectified  on-site, 
i.e.,  must  be  removed  to  MS  (here  assumed  to  occur  under  own 
propulsion),  but  at  (slow)  velocity  Vr  «  Vr-  As  will  be  discussed  later,  Rm 
might  include  repair  or  minor  failures  that  accumulated  during  the 
reverse  motion. 

{Note:  We  do  not  model  here  the  option  of  a  slow  progress  to  the  B 
for  repair). 

Rm(T(0))  =  Rm(N(T(0))):  Random  time  to  repair  all  minor  failures  (N(T(0)))  in 
number)  accumulated  during  generic  random  transit  time  T(0).  This  can 
represent  (a)  failures  of  sensors  (FLIR)  or  communications  equipment 
until  breakdown  (T(0)=Tb),  or  (b)  (any  such  “minor”)  failures 
accumulated  during  major-failure-fi'ee  transit  from  the  MS  to  the  B,  in 
time  t=d/VT. 

5.2  The  Basic  Stochastic  Model  for  Time  from  the  MS  to  the  B 

A  platform/vehicle  plans  to  travel  at  velocity,  Vj  from  the  MS  to  the  B,  a  distance  d. 
If  the  vehicle  breaks  down  before  covering  d,  it  proceeds  at  velocity  (Vr  «  Vr)  back  to 
the  MS,  where  it  is  repaired  and  re-launched. 
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Hence 


=  1  +  - 


s  +  R,,  +R 
M  m 


with  probability  Fj,  (ds),  0  <  s  <  t 
b 


where  we  assume  that  minor  failures  occur  Poisson-wise,  until  either  trip  completion  or  a 
serious  breakdown  occurs,  inducing  a  re-start.  Note:  It  has  been  tacitly  assumed  here  that 
minor  failure  repairs  at  the  MS  take  the  same  amount  of  time  as  at  the  B.  This  certainly 
need  not  be  true:  one  can  specify  Rm,  Ms(t)  and  R^,  b(s)  as  being  quite  different 
random  processes. 

In  (5.1),  Td'  refers  to  an  independent  random  replica  of  Td-  If  failure-inducing  fault 
removals  occur,  we  obtain  a  system  of  such  equations;  this  important  feature  is  omitted 
for  the  present. 

Also  in  (5.1),  the  accumulated  minor  failure  repair  times  can  be  modeled  as  a  random 
sum  of  independent  random  variables,  where  the  random  number  of  summands  depends 

upon  the  time  of  exposure  until  either  the  MS-B  transit  completes,  or  a  breakdown  and 
reverse  transit  occurs,  starting  at  2),.  For  time  of  exposure  T(0)=x,  and  {A^(j),  0<  y  <  x} 

a  Poisson  process  with  rate  and  Sj  a  random  minor  failure  service  duration, 

N(x) 

(5.2) 

J=0 

Comments:  The  above  model  formulation  is  illustrative  and  can  readily  be  changed, 
although  not  without  a  price  in  additional  model  complexity.  A  simpler  version  would 
ignore  all  failures  except  the  major  ones,  but  require  auxiliary  transport.  Such  a  version 


already  shows  the  amplifying  impact  of  in-transit  failures  on  overall  measures  of  transit 
time,  such  as  its  mean,  E\T^Y  This  effect  becomes  more  pronounced  if  it  becomes 

necessary  to  repair  accumulated  minor  failures.  It  has  also — and  importantly — been 
assumed  that  such  repairs  can  be  accomplished  flawlessly  during  the  visit  to  the  MS  to 
repair  any  major  failure,  but  clearly  this  may  not  occur.  Also,  in  this  model  no  attempt  is 
made  to  explicitly  represent  delays  in  queue  for  repair  at  the  MS.  These  delays  can 
introduce  a  substantial  increase  in  the  time,  Td,  and  hence  an  increase  in  the  time  imtil  a 
platoon  of  n  (here  treated  as  independent  and  identical  units)  assembles  at  the  B.  The 
reader  is  reminded  that  the  stochastic  behavior  of  the  individual  vehicles  (AVs)  in  a 
platoon  is  modeled  as  independent,  which  is  the  consequence  of  having  plentiful  support 
by  LCU(R)s  and  by  maintenance  facilities  at  the  MS  and  the  B.  Hence,  results  from  the 
present  model  are  almost  certainly  optimistic.  They  can,  however,  provide  useful  checks 
on  more  detailed  simulation  models. 

5.3  Expectation/Mean  of  Tj 

A  general  formula  for  the  expectation/mean  of  the  time  for  an  AV  to  complete  the 
transit,  T^,  is  obtained  by  conditional  expectation  from  (5.1);  first  condition  on  TJ,: 
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Consequently,  removing  the  condition  on  Tb, 


where 


q = [' + ^K(')](i  -  %  (')] + e[r^^]fj.^  (!) 


+  UFj-As) 


Cl  ='-%{<) 


5.5  Special  (Tractable)  Case:  Poisson  Minor  Failures 

For  illustration,  let  minor  failures  occur  according  to  a  stationary  Poisson  process, 
with  rate  Am,  with  repair  durations  {5^,7  =  1,2,... A''(7'(0))} independently  and  identically 

distributed,  and  with  mean  £^[*5]. 

The  result  is 


Note  that  the  expectations  with  respect  to  are  over  the  set  0  <  <  t ,  e.g.. 


E[Tb\t]  =  \‘xdFT,{xy, 
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this  applies  to  the  later  expectations,  £[^(7^);/],  as  well,  with  g(»)  being  an  arbitrary 
function  (of  T/,). 

5.6  Exponentially  Distributed  Tt, 

If  is  distributed  exponentially  with  rate  parameter  X  (a  very  special  case)  then 

E[Tb-,t]  =  \‘^xe-^^^Xdx  =  ^[l  -  (l  + 


and 


(5.7) 


l-Fn{t)  =  e-^. 


For  very  small  X  this  gives 


E[T,]  =  r[l  +  X„E[S]+XE[Rm^  +  (^1  +  ^ j{l  +  Xr„E[S'^ 


(5.8) 


AVs  in  transit  from  the  MS  to  the  B  may,  in  principle,  break  down  and  require 
transport  back  to  the  MS  for  repairs  several  times  before  successful  transit  is 
accomplished.  The  effect  on  eventual  transit  time  Td,  of  such  retum-to-MS  transports  and 
re-repairs  can  become  quite  large;  the  time  to  (slowly)  return  for  repairs  adds  to  the  total 
time,  and  the  repairs  congest  the  repair  facility,  adding  delay.  The  present  infant-failure 
model  does  not  permit  an  early  (infant)  failure  to  get  far  from  the  MS,  but  it  does 
compound  the  repair  burden  and  increases  turnaround  delay,  and  therefore  extends  the 
transit  time  to  the  B. 
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6.  Comparison  of  Deterministic  and  Stochastic  Models 

In  this  section  we  compare  results  from  the  stochastic  model  and  a  simplified 
deterministic  model.  The  deterministic  model  considered  in  this  section  is: 


^  dl  ^ 

Rate  of  "increase" 
of  number  of  active 
AVs  in  the  Littoral 


Rate  of  failure  Rate  of  "free"  active 
of  active  AVs  AVs  reach  ing  B 


(6.1) 


+ 


Ms(t) 


Rate  at  which  AVs  under 
maintenance  at  MS  become 
active  in  the  Littoral 


dQxCt) 

dl 

Rate  of  "increase"  of 
number  of  quiescent 
AVs  awaiting  tow 


TA(t) 

Rate  of  towable- 
faiiure  events 


^plQtCO 

V - ^ - / 

Rate  of  removal  initiation 
by  available  auxiliary 
vehicles  (e.g.,  LCs) 


(6.2) 


dRL(t) 

Rate  of  "increase"  of 
the  number  of  AVs 
being  transported  to 
the  MS  by  an  auxiliary 
vehicle 


VplQtW 

Rate  of  pickup  of  failed 
AVs  by  available  landing 
craft;  headed  for  the  MS 


Vri^RlW 

V . .  . 

Rate  of  delivery 
of  landing  craft 
towed  AVs  to 
the  MS 


(6.3) 


dM3(t) 

Rate  of  "increase" 
of  number  of  AVs 
in  MS  maintenance 


V'rlRl(I) 

V - ^ - / 

Rate  of  dropoff  at 
MS  from  LC 
transport 


1  +  Ms(t) 

V - V  ' 

Rate  of  maintenance 
completion  at  the  MS 


(6.4) 


dB(t) 

Rate  of  "increase" 
of  active,  available 
AVs  at  the  B 


=  vvA^ 

Rate  of  "free’ 
active  AVs 
reaching  the 


(6.5) 
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In  this  model,  all  disabled  AVs  are  towed  back  to  the  MS.  There  is  always  an  LC 
available  for  towing. 

The  corresponding  stochastic  model  for  one  AV  is  that  of  (5.1),  with  the  minor  failure 
rates  Am  =  0.  The  variance  of  the  time  for  a  single  AV  to  arrive  at  the  B  for  the  stochastic 
model  of  this  section  can  be  computed  using 


[{dlVjf 


if  T^>d/VT 


Vr  Ft' 

y 

Tb  +  Rm  +  Td 

.IL 

J 

if  n<d/VT 


A 


where 


(6.6) 


1  + 


Vj 

v'r\ 


2  dlVj 

I  x^FT,{dx)+2[E[Rm]^E[TA 

0  I 


xFn{dx) 


+[e[rI]  +  2E[Rm]E[Td'^lFT, 

andFr^(r)  =  l-Fr*(t) 


We  approximate  the  distribution  of  the  AV  arrival  time  at  the  B  for  an  exponential 
time  to  failure  as  follows:  if  no  failure  occurs,  then  the  time  for  an  AV  to  reach  the  B  is 
dWt  hours;  the  probability  of  at  least  one  serious  failure  in  transit  is 
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Then  we  approximate  the  distribution  of  the  time  for  an  AV  to  reach  the  B  as 


X  =  —+y  where 
Vt 


Z  withprob/7/ 

0  with prob  I- Pf 


and  Z  has  an  exponential  distribution  with  mean 


E[Z]  = 


Vt) 


ipf- 


Note  that 


The  variance  for  the  approximation  (6.7)  is 


(6.7) 


(6.8) 


Var[X]  =  Var[Y]  =  Var[Z]pf  +  E[zf  Pf{l-pf) 


(6.9) 


The  stochastic  model  corresponding  to  the  deterministic  model  for  multiple  vehicles 
is  that  each  of  the  vehicles  travels  independently  to  the  B;  all  the  AVs  start  at  time  zero. 
Thus  the  number  of  vehicles  that  have  arrived  at  the  B  by  time  s  has  a  binomial 
distribution,  with  number  of  trials  being  the  number  of  AVs,  and  the  probability  of 
success  being  P{X  <  5}. 

6.1  Numerical  Example 

Below  are  displayed  results  for  a  numerical  example.  The  parameter  values  are: 

Deterministic  Model: 

Rate  of  transit  completion  to  the  beach  for  a  working  AV  is  Vf  =  20  /  25  per  hour 

Rate  of  commencing  towing  for  failed  AV  is  vpi  =2  per  hour 

Rate  of  tow  completion  to  mother  ship  (MS)  is  Vrl  =  1/5  per  hour 

Rate  of  repair  for  failed  vehicle  at  MS  is  ps  =  0.5  per  hour 

Rate  of  failure  for  a  working  AV  is  /I  =  1  /  3  per  hour 
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Stochastic  Model: 

Distance  to  B,  d  is  25nm 

Planing  velocity.  Ft,  for  a  working  AV  is  20nm  per  hour 
Towing  velocity,  F/e,  to  the  MS  is  5nm  per  hour 
Mean  repair  time  at  the  MS  is  two  hours 
Variance  of  a  repair  time  at  the  MS  is  one  hour^ 

Mean  time  to  serious  failure  is  three  hours 


We  assume  first  that  the  time  to  a  serious  failure  has  an  exponential  distribution.  If  no 
serious  failure  occurs,  then  the  time  for  an  AV  to  reach  the  B  is  1.25  hours.  In  the 
stochastic  model,  allowing  for  failures  as  above,  the  mean  time  for  an  AV  to  reach  the 
beach,  Tj,  is  3.79  hours,  and  the  standard  deviation  of  the  time  is  4.59  hours.  The 
probability  of  at  least  one  serious  failure  in  transit  is  pf  -  I-  expj— jl.25|  =  0.34 . 

We  will  approximate  the  distribution  of  the  time  for  an  AV  to  reach  the  B  as 


X  =  125  +  Y  where 

iZ  with  prob  p^ 

0  with  prob  1-Pf 


(6.10) 


and  Z  has  an  exponential  distribution  with  mean 


E[Z]  =  {E[Td]-U5)/pf.  (6.11) 

The  mean  of  the  approximate  time  is  3.79  hours,  the  same  as  that  of  Tj.  However,  its 
standard  deviation  is  5.60  hours,  which  is  22%  greater  than  that  of  Td. 
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Figure  4  below  displays  the  mean  and  quantiles  for  the  binomial  stochastic  model  and 


the  mean  from  the  deterministic  model. 


Time  (Hours  from  MS  Launch) 

r -  Fluid  Modd  ~  ~ 

;  .  .  .  Binomial  Model  (mean) 


- Binomial  Model:  smallest  k  s.t.  Prob(B(t)<=k)  >=0.05 

- Binomial  Model:  smallest  k  s.t.  Prob(B(t)<=k)  >=0.95 


Figure  4. 

Number  of  AVs  at  the  Beach  Time  from  Launch. 


Comments:  It  can  be  seen  from  Figure  4  that  the  deterministic  and  stochastic  model 
means  agree  well  when  the  time  exceeds  one  hour.  The  stochastic  model  permits 
probability  limits  to  be  created;  these  suggest  considerable  variability  in  the  number  of 
platoon  elements  (AVs)  that  have  arrived  at  the  B  by  any  fixed  time  t.  For  instance,  at 
f  =  6  hours,  the  probability  that  eight  or  more  out  of  12  have  arrived  is  ~  95% ,  while  the 
mean  number  is  =  10.  It  is  likely  that  such  variability  is  actually  too  small. 
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In  the  results  displayed  below  we  explore  the  effect  of  different  distributions  of  the 
time  to  failure  on  the  time  for  an  AV  to  arrive  at  the  B.  We  use  the  stochastic  model  of 
this  section.  We  compare  the  mean  and  standard  deviation  of  the  arrival  time  at  the  B,  Tj, 


for  times  to  serious  failure  having  a  mixed  exponential  distribution:  for  x  >  0, 

The  following  tables  display  the  mean  and  standard  deviations  for  the  AV  arrival 
times  at  the  B.  Note:  Once  again,  the  form  of  the  distribution  of  the  time  to  serious 
failure  can  have  a  dramatic  effect  on  the  mean  and  standard  deviation  of  the  time  to  arrive 
at  the  B.  The  form  of  that  distribution — specifically  its  tendency  to  exhibit  small 
values — will  also  influence  logistics  requirements,  and  their  spatial  location.  Comparison 
of  the  mean  and  standard  deviation  of  the  operationally  important  time  to  arrive  at  the  B, 
suggests  that  the  eventual  arrival  time  distribution,  Td,  will  be  approximately  exponential 
in  the  case  of  frequent,  serious  failures  of  the  AV. 
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Distribution  of  Time 
to  Failure 

Velocity 
When  not 
Failed 

Time  to 
Arrive  at 
Destination 
if  No  Failure 

Mean  Time 
to  Arrive  at 
Destination 

Standard 
Deviation 
of  Time  to 
Arrive  at 
Destination 

(Hours) 

(Hours) 

(Hours) 

Exponential  Mean  3 

5 

5.00 

29.36 

27.69 

10 

2.50 

9.31 

9.45 

15 

1.67 

5.40 

mmm 

1.25 

3.79 

4.59 

25 

1.00 

2.91 

mmm 

30 

0.83 

2.36 

3.30 

35 

1.99 

2.94 

40 

1.72 

2.68 

Mixed  Exponential  Mean  3 
p=0.5;  1,  X2~1/5 

5 

5.00 

25.85 

23.81 

10 

2.50 

15 

1.67 

7.91 

20 

1.25 

5.29 

6.21 

25 

1.00 

30 

0.83 

3.41 

4.53 

35 

0.71 

2.89 

40 

0.63 

2.50 

3.68 

Table  1. 

Moments  for  Time  of  Individual  AV  to  Arrive  at  Destination,  Beach  (B),  which  is 
25nm  away.  The  velocity  of  return  to  Mother  Ship  (MS)  if  failed  is  five  nm  per 
hour;  the  mean  repair  time  is  two  hours  and  the  variance  of  the  repair  time  is 
one  hour. 
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Distribution  of  Time 
to  Failure 

Velocity 
When  not 
Failed 

Time  to 
Arrive  at 
Destination 
if  No  Failure 

Mean  Time 
to  Arrive  at 
Destination 

Standard 
Deviation 
of  Time  to 
Arrive  at 
Destination 

(nm/Hour) 

(Hours) 

(Hours) 

(Hours) 

Exponential  Mean  6 

5 

5.00 

13.21 

11.45 

10 

2.50 

5.34 

5.15 

15 

3.33 

3.58 

20 

2.41 

2.86 

25 

1.00 

2.43 

30 

0.83 

2.14 

35 

0.71 

1.32 

1.94 

1.15 

1.78 

Mixed  exponential  Mean  6 
p=0.5;?ii=l,  >.2=1/11 

5 

||||■|■||^ 

I^B 

n 

2.50 

8.32 

8.17 

15 

1.67 

5.81 

6.33 

20 

1.25 

4.45 

5.24 

25 

1.00 

3.56 

4.51 

30 

0.83 

2.97 

4.00 

35 

2.54 

3.63 

40 

0.63 

2.22 

3.33 

Table  2. 

Moments  for  Time  of  Individual  AV  to  Arrive  at  Destination,  Beach  (B),  which  is 
25nm  away.  The  velocity  of  return  to  Mother  Ship  (MS)  if  failed  is  five  nm  per 
hour;  the  mean  repair  time  is  two  hours  and  the  variance  of  the  repair  time  is 
one  hour. 
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Distribution  of  Time 
to  Failure 

Velocity 
When  not 
Failed 

Time  to 
Arrive  at 
Destination 
if  No  Failure 

Mean  Time 
to  Arrive  at 
Destination 

Standard 
Deviation 
of  Time  to 
Arrive  at 
Destination 

(nm/Hour) 

(Hours) 

(Hours) 

(Hours) 

Exponential  Mean  9 

5 

5.00 

9.86 

7.89 

10 

2.50 

4.29 

3.87 

15 

1.67 

2.73 

2.77 

20 

1.25 

2.00 

2.24 

25 

1.00 

1.58 

1.92 

30 

0.83 

1.31 

1.70 

35 

0.71 

1.11 

1.54 

40 

0.63 

0.97 

1.42 

Mixed  exponential  Mean  9 
p=0.5;  Xi=l,  ^2=1/17 

5 

5.00 

12.48 

10.04 

10 

2.50 

7.52 

7.30 

15 

1.67 

5.41 

5.88 

20 

1.25 

4.18 

4.95 

25 

1.00 

3.39 

4.32 

30 

0.83 

2.84 

3.85 

35 

0.71 

2.44 

3.50 

40 

0.63 

2.14 

3.22 

Table  3. 

Moments  for  Time  of  Individual  AV  to  Arrive  at  Destination,  Beach  (B),  which  is 
25nm  away.  The  velocity  of  return  to  Mother  Ship  (MS)  if  failed  is  five  nm  per 
hour;  the  mean  repair  time  is  two  hours  and  the  variance  of  the  repair  time  is 
one  hour. 
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Distribution  of  Time 
to  Failure 

Velocity 
When  not 
Failed 

Time  to 
Arrive  at 
Destination 
if  No  Failure 

Mean  Time 
to  Arrive  at 
Destination 

Standard 
Deviation 
of  Time  to 
Arrive  at 
Destination 

(nm/Hour) 

(Hours) 

(Hours) 

(Hours) 

Exponential  Mean  1 8 

5 

7.17 

4.73 

10 

3.34 

2.52 

15 

1.67 

2.18 

1.85 

20 

1.25 

1.62 

1.51 

_ 25 

1.28 

1.17 

35 

0.71 

0.63 

0.79 

0.98 

Mixed  exponential  Mean  1 8 
p=0.5;Xi=l/30,  >12=1/6 

5 

5.00 

8.61 

6.43 

10 

2.50 

3.97 

3.44 

15 

1.67 

2.57 

2.52 

20 

1.25 

1.90 

2.06 

25 

1.51 

1.78 

30 

1.25 

35 

1.06 

lllllllllllll^^ 

40 

0.63 

0.93 

1.33 

Table  4. 

Moments  for  Time  of  Individual  AV  to  Arrive  at  Destination,  Beach  (B),  which  is 
25nm  away.  The  velocity  of  return  to  Mother  Ship  (MS)  if  failed  is  five  nautical 
miles  (nm)  per  hour;  the  mean  repair  time  is  two  hours  and  the  variance  of  the 
repair  time  is  one  hour. 
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